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Abstract 

Plasmodium falciparum malaria imposes a serious public health concern throughout the tropics. Although 
genetic tools are principally important to fully investigate malaria parasites, currently available forward and 
reverse tools are fairly limited. It is expected that parasites with a high mutation rate can readily acquire novel 
phenotypes/traits; however, they remain an untapped tool for malaria biology. Here, we generated a mutator 
malaria parasite (hereinafter called a 'malaria mutator'), using site-directed mutagenesis and gene transfec- 
tion techniques. A mutator Plasmodium berghei line with a defective proofreading 3' ■ 5' exonuclease activ- 
ity in DNA polymerase 8 (referred to as PbMut) and a control P. berghei line with wild-type DNA polymerase 5 
(referred to as PbCtI) were maintained by weekly passage in ddY mice for 122 weeks. High-throughput 
genome sequencing analysis revealed that two PbMut lines had 1 75-1 78 mutations and a 86- to 90-fold 
higher mutation rate than that of a PbCtI line. PbMut, PbCtI, and their parent strain, PbWT, showed similar 
course of infection. Interestingly, PbMut lost the ability to form gametocytes during serial passages. We 
believe that the malaria mutator system could provide a novel and useful tool to investigate malaria biology. 
Keywords: mutator; Plasmodium; DNA polymerase 8; genome sequencing 
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deaths worldwide in 201 0, due mostly to infections of 
Plasmodium falciparum, the most virulent human 
malaria parasite. 1 This parasite has developed drug re- 
sistance to almost all conventional antimalarial drugs. 2 
Development of malaria vaccines targeting both insect 
and blood stage parasites is still ongoing. 3,4 These 
circumstances urge us to further pursue effective in- 
tervention strategies to combat malaria. However, our 
knowledge of the parasite biology to identify new 
drug and vaccine targets is still inadequate. Current 
genetic tools, although limited, provided opportunities 
to specifically alter the parasite genome to explore its 
biology and gain new insights into gene function and ex- 
pression. 5-1 1 A combination approach that exploitsthe 
advantages of different genetic tools may provide the 
best strategy to study malaria parasite biology. 

Important biological phenotypes/traits of malaria 
parasites like drug resistance are often caused by a 
point mutation(s). 12-14 If the mutation rate of 
malaria parasites is artificially increased, it is expected 
that parasites can rapidly acquire novel phenotypes/ 
traits that can be selected for under appropriate select- 
able conditions: e.g. drug resistance would be readily 
produced under the pressure of a drug of interest. To 
date, little attempts have been made to utilize such 
parasites with an accelerated mutation rate to study 
malaria parasite biology. 

In eukaryotic organisms, the spontaneous mutation 
rate is very low, e.g. 3.3 x 1 O~ 10 /basepair(bp)/celldiv- 
ision in yeast, 15 2.7 x 1 0~ 9 /bp/generation in Caenor- 
habditis eleqans^ 6 and 5.8 x 1 0~ 9 /bp/generation in 
Drosophila meianoqaster} 7 The mutation rate of in 
vitro P. falciparum culture was recently reported as 
1.0-9.7 x 1 0~ 9 /bp/generation. 1 8 These low muta- 
tion rates are ensured by the high fidelity of DNA repli- 
cation due to the strict selectivity of correct bases and 
the proofreading function of DNA polymerases, as 
well as the mismatch repair system. DNA polymerases 
a, 8, and e are involved in nuclear genome DNA synthe- 
sis, ofwhich polymerases 8 and e possessa proofreading 
3' — > 5' exonuclease activity for lagging and leading 
DNA strand replication, respectively. The proofreading 
activity domain of DNA polymerases 8 and e is widely 
conserved from prokaryotestoeukaryotes, 1 9,20 includ- 
ing Plasmodium. 2 ^ Disruption of the proofreading activ- 
ity by site-directed mutagenesis significantly increased 
a mutation rate by as much as 10-fold to >1 00- 
fold, 20,22,23 and the mutator effect is particularly 
strong in DNA polymerase 8. 20,24 Furthermore, in 
microorganisms such as yeast and fungi, the mutator 
system has been employed to generate various mutated 
phenotypes including drug resistance. 25-28 

In the present study, we generated a mutator malaria 
parasite (hereinafter called a 'malaria mutator') using 
a rodent malaria parasite, Plasmodium berghei ANKA 
by eliminating the proofreading activity of its DNA 



polymerase 8 and performed genome-wide analysis of 
mutations in order to see whether mutation rate in 
the mutator increased as expected. The mutation rate 
of the malaria mutator was 86- to 90-fold higher 
than the control parasite. 

2. Materials and methods 

2.1 . Generation of mutator parasites and its phenotypic 
analyses 

The rodent malaria parasite, P. berghei ANKA strain 
(clone 2.34), was used. The plasmid and primer 
sequences used for generation of mutator parasite are 
described in Supplementary Data (Supplementary 
TableS! ). In brief,thetwoaminoacidscriticalforproof- 
reading in DNA polymerase 8 were replaced with 
alanine in the mutator parasites. The phenotype ana- 
lyses (growth rate, gametocytogenesis, male gameto- 
genesis, and ookinete formation) of the mutator 
parasites were performed as described in Supplemen- 
tary Data. 

2.2. Maintenance of parasites in mice 

Mutator (referred to PbMut) and control parasites 
(PbCtl) were maintained by weekly passages in mice. 
We injected intraperitoneal^ 100 or 1,000 infected 
erythrocytes into two female ddY mice (6- to 10- 
week-old; Japan SLC), respectively. Infected blood of 
either mouse reaching adequate parasitaemia (several 
percent) was collected 1 week later by cardiac puncture 
and used for the next passage. Cloning of parasites by 
limiting dilution was done in the 53rd and 119th 
week of passage for PbMut, and in the 43rd and 
1 1 9th week of passage for PbCtl. 

2.3. High-throughput genome sequencing 
Genomic DNA from the cloned lines and PbWT 

(Supplementary Data) was sequenced genome wide 
using an lllumina HiSeq 2000 or an lllumina Genome 
Analyzer II. Paired-end read sequencing was performed 
usingread lengthsof 90or 1 00 bases. Lowquality reads 
with an N-rate of >1 0% and/or reads for >50% all 
bases with a PHRED quality score of <5 were removed. 
To exclude mouse DNA-derived reads, the sequence 
reads were aligned to the Mus musculus (C57BL/6J) 
genome sequence from the NCBI Genome database 
using SOAP2, 29 and matched reads were removed. 

2.4. Read mappinq and mutation identification 

The sequence data were further processed using two 
pipelines, i.e. one by read alignment using Bowtie 
version 0.1 2. 7 30 and variant calling with VarScan 
version 2.2.8, 31 and the other with the CLC Genomics 
Workbench software (CLC Bio). We used the genome 
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sequence (1 8.2 Mb) of the P. berghei AN KA strain 32 as a 
reference sequence, which was downloaded from the 
Wellcome Trust Sanger Institute website (http://www. 
sanger.ac.uk/). Further details are provided in Supple- 
mentary Data. Mutations detected by lllumina genome 
sequencing were validated preliminarily by Sanger se- 
quencing. For comparison, all insertions and deletions 
(indels), all base substitutions causing non-sense muta- 
tions,and all base substitutions detected usingeitherof 
two pipelines were validated preliminarily by Sanger 
sequencing (Supplementary Table S2). The Sanger se- 
quencing confirmed all the mutations (data not 
shown). Thus, we inferred thatthe sequences produced 
by lllumina sequencing were consistent. 

2.5. Detection of mutation clusters in a chromosome 
The r-scan statistics was used to detect significant 

clustering of mutations along each chromosome as 
described previously. 33 We used two probability cut- 
offs: 1 and 5%. Overlappingclustering regions were con- 
sidered as a single mutation cluster region. 

2.6. Estimation of mutation rates 

We estimated the mutation rate using the base sub- 
stitution data generated by the Bowtie and VarScan 
pipeline. The mutation rate was calculated using the 
equation jjL bs = m/(«T), where |x bs is the base substitu- 
tion rate (per nucleotide site per day), m is the 
number of mutations detected, n is the number of nu- 
cleotide sites, and T corresponds to the passage 
days. 34 The mutation rate was ca Iculated for the follow- 
ing specific regions: synonymous sites (Syn) and non- 
synonymous sites (NonSyn) in protein-coding regions; 
introns; intergenic regions; and the overall sequences. 
Syn and NonSyn sites were estimated as described by 
Nei and Gojobori. 35 The mutation rate per replication 
was also calculated. In the ddY mice used here, P. 
berghei multiplied ~1 0-fold perdayduringthe logarith- 
mic growth phase, and hence, we assumed log 2 10 
(«3.3) replications per day. We considered that the 
number of base substitutions were the same as that of 
point mutations because no complicated mutations 
were found at a given substituted site. Fisher's exact test 
for the 2x2 contingency table was performed to 
assess statistical significance of difference in variant pro- 
portion between specific regions (i.e. Syn vs. NonSyn, 
NonSyn vs. intron, and intron vs. Syn). For this test, we 
rounded the numbers of Syn and NonSyn sites to unit. 

2.7. Data access 

The raw sequence reads used in this study are available 
from the DDBJ Sequence Read Archive (DRA) under 
accession No. DDBJ: DRA000656. The information for 
P. berghei genes discussed in this publication is available 
from PlasmoDB (www.plasmodb.org). 



3. Results 

3.1 . Generation of malaria mutator 

Using site-directed mutagenesis and gene transfec- 
tion techniques, we established a mutator P. berghei 
line where DNA polymerase 8 had a defective proof- 
reading activity (referred to as PbMut) and a control P. 
berghei line with wild-type (non-mutated) DNA poly- 
merase 8 (referred to as PbCtl) (Supplementary Fig. 
S1 ). These parasite lines were then subjected to serial 
passages and cloning, and two PbMut lines (M1 22A 
and M1 22B) and one PbCtl line (C1 22), as well as one 
PbWT line, were obtained (Supplementary Fig. S2). 

3.2. Mutations generated in the malaria mutator 
High-throughput genome sequencing of PbMut lines 

(M1 22Aand M1 22B),one PbCtl line (C1 22), and one 
PbWT line produced each 0.9-1.2 Gb, which corre- 
sponded to ~ 50-fold redundancy relative to the P. 
berghei ANKA genome sequence from the Wellcome 
Trust Sanger Institute (1 8.2 Mb) (Supplementary 
Table S3). We reconstructed ~86-95% sequence 
regions that encompassed the 14 chromosomes and 
analysed them for mutations generated in PbMut and 
PbCtl. Preliminary analysis revealed that PbWT (clone 
2.34) used in this study had 71 substitutionscompared 
with the Wellcome Trust Sanger ANKA reference se- 
quence (Supplementary Table S4). It is not certain 
whether these differences are true substitutions or 
due to artefacts such as ma pping errors and/or sequen- 
cing/assemblyerrorsfrom the lowercoverage (x4) ref- 
erence genome sequence. 32 We therefore compared 
the sequence of PbWT (clone 2.34) with those of 
PbMut and PbCtl lines, and different sites we re regarded 
as variants that occurred in each line. The number of 
point mutations in PbMut and PbCtl differed dramatic- 
ally (Table 1, Fig. 1, and Supplementary Table S5). 
We only detected two base substitutions in C122, 
whereas 1 78 and 1 75 basesubstitutionsweredetected 
in M122A and M122B, respectively (Supplementary 
Table S5). A large proportion of the substitutions (n = 
1 62) were shared by the two PbMut lines. This high 
rate of shared mutations between the two PbMut 
clones was probably due to a severe bottleneck effect 
during each mouse passage, where 100-1,000 
infected erythrocytes were inoculated intraperitone- 
al^, corresponding to intravenous injection of 1 0- 
1 00 infected erythrocytes. There were 1 07 and 1 1 0 
base substitutions in the protein-coding gene regions 
of M1 2 2A and M1 2 2B, respectively. Most of these sub- 
stitutions, 89 and 92 substitutions in M1 22A and 
M1 22B, respectively, were NonSyn substitutions (mis- 
sense mutations). 

There were five non-sense mutations: in the serine/ 
threonine protein kinase gene (PBANKA_090380), the 
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Table 1. Summary of substitutions in the P. berghei control and mutator clones 3 





Control 




Mutator 






C1 22 b 




M1 22A b 


M1 22B b 


Total number of sites 


1 6.5M 




16.1M 


1 6.6M 


Number of base substitutions 


2 




1 76 


1 73 


Total number of sites in protein-coding regions 


9.7M 




9.6M 


9.7M 


Number of substitutions in protein-coding regions 


2 




1 05 


108 


dS c 


0 




1.1 x 1 cr 8 


1.1 x 10~ 8 


Total number of Syn sites 


1.9M 




1.9M 


1.9M 


Number of Syn substitutions 


0 




1 8 


1 8 


dN c 


2.9 x 1 (T 10 




1.3 x 1 cr 8 


1.3 x 10~ 8 


Total number of NonSyn sites 


7.8M 




7.7M 


7.8M 


Number of NonSyn substitutions 


2(0) d 




87(3) d 


90(3) d 


Total number of sites in introns 


1.1M 




1.1M 


1.1M 


Number of substitutions in introns 


0 




8 


7 


Total number of sites in intergenic regions 


5.6M 




5.4M 


5.7M 


Number of substitutions in intergenic regions 


0 




63 


58 


a Base substitutions detected from alignments by the Bowtie software are listed. 

b Numbers followed by 'C (control) and 'M' (mutator) indicate the number of weeks of mouse passaging for each clone. 

c dS and dN are the number of Syn substitutions per Syn sites and the number of NonSyn substitutions per NonSyn sites, respectively. 

d ln parenthesis, the number of substitutions resulting in non-sense mutation is shown. 
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Figure 1. Distribution of mutations on the 1 4 chromosomes of P. berghei control and mutator lines. (A) Mutations in a P. berghei control line 
(C1 22). (B) Mutations in two P berghei mutator lines (the upper and lower horizontal bars represent lines M1 22Aand M1 22B, respectively). 
The 1 4 chromosomes are represented by horizontal bars with grey and dark blue blocks. The grey blocks indicate the sequence regions used 
for mapping to detect mutations, and the dark blue blocks were those regions not used for mapping because of their lower coverage (see 
Materials and Methods for details). The red vertical bars indicate point mutations detected using the Bowtie program. The cyan vertical 
bars indicate short indels detected by the CLC Genomics Workbench software. The green and yellow blocks indicate clusters of 
mutations with the significance levels of P < 0.01 and <0.05, respectively. Asterisks show closely spaced vertical bars that could not be 
separated at the resolution scale used for this figure. 



secreted ookinete protein gene (PBANKA_1 23360), in introns and 58-63 base substitutions in intergenic 
and three genes with unknown functions (Table 2). regions. PbCtl had two NonSyn substitutions in genes 
PbMut lines also had seven to eight base substitutions with unknown function (Supplementary Table S5). Five 
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Gene ID 


Product (total amino acid length) 


AA position 


Clone 


PBANKA_083220 


Conserved Plasmodium protein, unknown function (300) 


E87 


M1 22B 


PBANKA_090380 


Serine/threonine protein kinase (965) 


E677 


M1 22A a ,M1 22B a 


PBANKA_1 1 3490 


Conserved Plasmodium protein, unknown function (870) 


E336 


M1 22A.M1 22B 


PBANKA_1 23360 


Secreted ookinete protein (205) 


E144 


M1 22A 


PBANKA_1 23 730 


Conserved Plasmodium protein, unknown function (486) 


S465 


M122A, M122B 



a Base substitutions were detected using CLC Genomics Workbench. 

Table 3. Indel mutations in P. berghei mutator clones 



Chr Position Change 



Region 



Gene ID 



Product 



Clone 



4 


446,1 69 


(A) 4 - 


>(A) 3 


Interj 


»enic 






M1 22A,M1 22B 


8 


941 ,899 


(T) 10 - 


-(T)„ 


Interg 


»enic 






M1 22A,M1 22B 


1 3 


30,91 7 


(A) 2 - 


>(A), 


Interj 


»enic 






M122A, M122B 


1 3 


2,414,1 37 


COs- 


(T) 2 


Exon 




PBANKA_1 36380 


Phosphatidylinositol transfer protein 


M122A, M122B 


14 


1,372,440 


(c) 3 - 


>(.C) 2 


Exon 




PBANKA_1 43750 


Transcription factor with AP2 domain(s) (ApiAP2) 


M1 22A.M1 22B 



Indels were detected using CLC Genomics Workbench. 



single-base indel mutations were detected in single- 
nucleotide repeat regions in PbMut lines (Table 3), but 
none in the PbCtl line. Two were found in protein- 
coding genes, causing frameshifts in the phosphatidyl- 
inositol transfer protein gene (PBANKA_1 36380) and 
a gene encoding for ApiAP2 (PBANBKA_1 43 750). 

Analysis of chromosomal distribution of the muta- 
tions by r-scan statistics 33 revealed that there were 
nine chromosomal regions, where mutations were sig- 
nificantly clustered (Fig. 1) (P<0.05). This might 
suggest the presence of mutation hot spots in the P. 
berghei genome. 

3.3. Mutation rate of the mutator 

Overall substitution rate (mean + SE) was 0.14 x 
10" 9 /bp/day (0.42 x 1 0" 10 /bp/cell division) for 
C1 22, 1 2.7 x 1 0" 9 /bp/day (3.8 x 1 0~ 9 /bp/cell div- 
ision) for M122A, and 12.2 x 1 0" 9 /bp/day (3.7 x 
1 0~ 9 /bp/cell division) for M1 22B, respectively. Thus, 
the mutation rate was 86- to 90-fold higher in PbMut 
than in PbCtl. Since no mutation was detected in the 
intergenic regions in PbCtl, we were unable to calculate 
the mutation rate in the intergenic regions, but it 
should be <0.21 x 1 0" 9 /bp/day (0.62 x 1 0" 10 /bp/ 
cell division), the rateestimated byassumingone muta- 
tion inthe regions. The mutation rates (M1 22A: 1 3.6 x 
10" 9 /bp/day and M122B: 1 1 .9 x 1 0" 9 /bp/day) in 
the intergenic regions in PbMut were similar to those 
in overall regions (Fig. 2). In intragenic regions, the mu- 
tation rate in NonSyn sites was comparable with that of 
overall sites. Although the NonSyn mutation rates 
(M1 22A: 1 3.2 x 1 0" 9 /bp/day and M1 22B: 1 3.4 x 
1 0~ 9 /bp/day) appeared to be higher than the rates in 



X 

' CO 



B 




Syn NonSyn 



Intron 



IGR 



Overall 



Figure2. Mutation rate of P. berghei mutator lines. The mutation rate 
(per nucleotide site per day) was calculated for the following 
specific sequence regions: Syn and NonSyn sites in protein- 
coding regions; introns; IGR; and the overall sequences. A and B 
indicate Ml 22A and M1 22B, respectively. This figure appears in 
colour in the online version of DNA Research. 

Syn sites (M1 22A: 1 1.3 x 1 0" 9 /bp/day and M1 22B: 
1 1 .2 x 1 0" 9 /bp/day) and in introns (M1 22A: 8.5 x 
10" 9 /bp/day and M122B: 7.3 x 1 0" 9 /bp/day), 
these differences were not significant (P>0.05 by 
Fisher's exact test) (Fig. 2). 



3.4. Growth rate and sexual differentiation 
of the mutator 
The growth rate analysis revealed that two PbMut 
lines (M1 22A and B) proliferated in the mouse with a 
similar rate to those of PbCtl (C122) and PbWT 
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(Supplementary Fig. S3). This result suggests that elim- 
ination of the proofreading activity of DNA polymerase 
8 and the generated mutations did not affect asexual 
blood stage parasite growth. It seems likely that para- 
sites with mutations deleterious to parasite survival 
were eliminated during the propagation of parasites 
in mice, and thus parasites that may have showed slow 
growth rate were masked underthe presentexperimen- 
tal conditions. Interestingly, M122A and M122B did 
not produce gametocytes, while C122 did (Supple- 
mentary Table S5). Consistently, neither exflagellation 
nor ookinete formation was observed in PbMut lines 
(Supplementary Table S5). 



4. Discussion 

The present study demonstrated that elimination of 
the proofreading activity of the P. berghei DNA polymer- 
ase 8 produced malaria parasites (malaria mutator) 
with a high mutation rate (86- to 90-fold higher muta- 
tion rate in PbMut vs. control PbCtl).The increased mu- 
tation rate in PbMut is comparable with a >60 time 
increase observed in Saccharomyces cerevisiae. 23 We 
consider this 86- to 90-fold higher mutation rate as a 
tentative one that needs to be refined by further inves- 
tigation. In this study, only one parasite lineage each for 
PbCtl and PbMut was used to calculate the relative mu- 
tation rate. To obtain more accurate mutation rates in 
both PbCtl and PbMut, mutation rates from several 
lineages would be required. 

PbMut lines grew with a similar rate to those of PbCtl 
and PbWT lines (Supplementary Fig. S3). Because of 
this, it is unlikely that the 89-92 NonSyn mutations 
(and other mutations) in the PbMut lines were del- 
eterious to parasite survival. Apart from genes with 
unknown function, most of the genes with NonSyn 
changes code for housekeeping functions but not pro- 
teins involved in immune evasion (Supplementary 
Table S5). These genes may well have had relaxed con- 
straints, not seriously affecting parasite survival during 
the asexual blood stage. We identified five non-sense 
mutations in PbMut lines (Table 2),causingtruncations 
in each gene. The truncation of these genes, likewise, 
appeared not to affect parasite growth, suggesting 
that expression of these genes may not be essential for 
the asexual blood stage development. This is the case 
for the ookinete protein gene (PBANKA_1 23360), 
which is not expressed at the blood stage (PlasmoDB). 
Previously, it was shown that knockout of the serine/ 
threonine protein kinase gene (PBANKA_090380) was 
lethal, which suggested the importance of this gene 
during the asexual blood stage. 36 However, the mutation 
detected in PbMut lines was located downstream of the 
protein's kinase domain, thus the kinase activity was 
probably not impaired in the truncated protein. 
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Two indel mutations causing frameshifts were also 
detected: one in the phosphatidylinositol transfer 
protein gene (PBANKA_1 36380) and the other in the 
transcription factor gene with AP2 domain(s) (ApiAP2) 
(PBANKA_1 43750). The frameshift mutation at amino 
acid position 506 caused a truncation (1513 amino 
acids) in the Pleckstrin homology domain of the phospha- 
tidylinositol transfer protein gene. This domain binds to 
phosphatidylinositol lipids, suggesting that the truncated 
protein lost its function. BLASTP search did not identify 
any homologue of PBANKA_1 36380 in PlasmoDB. Thus, 
the expression of the phosphatidylinositol transfer 
protein gene may not be essential during the asexual 
blood stage. The ApiAP2 gene is a member of the gene 
family of DNA binding transcription factor proteins and 
is involved in the regulation of gene expression at various 
stages of parasite development. 37 The frameshift in the 
ApiAP2 gene caused an extension of the gene product 
(from 2339 to 2377 amino acids in non-mutated gene 
and frameshift- mutated gene, respectively). Recently, it 
was reported that two clones of gametocyte-deficient P. 
falciparum produced by piggyBac random mutagenesis 
were attributed to insertions in ApiAP2 gene. 38 Notably, 
we did detect no gametocyte in PbMut lines at the 
1 22nd mouse passage. The ApiAP2 gene which we identi- 
fied is not orthologous to one that Ikadai et al. reported, 
but the APiAP2 family may function as a key molecule 
for sexual differentiation of Plasmodium parasites. 

Random mutagenesis using artificial mutagenssuchas 
alkylating agents or irradiation remains viable tools to 
generate malaria mutants. However, artificial mutagen- 
induced mutationsare highly biased, e.g.G/CtoA/Ttran- 
sition by ethyl methanesulfonate, 25,39 A to T and T to A 
transversions by N-ethyl-N-nitrosourea, 39 or pyrimidine 
dinucleotide-associated C to T transitions by UV or 
neutron irradiation. 40 These biases may influence 
studies such as drug resistance and vaccine resistance, 
which should reproduce the actual mutants in the field. 
The mutations generated in the malaria mutator were 
probably generated during natural processes of DNA rep- 
lication as the principle of the mutator relies on a defect- 
ive proofreading activity of DNA polymerase 8. Malaria 
parasites with identified DNA mismatch repair activities 
have been linked to drug-resistant parasites. 41 Because 
of this, we believe that the malaria mutator system 
could be useful for reproducing drug resistance observed 
in the field. We have, in fact, generated mutator parasites 
resistant to sulfadoxine and could identify a NonSyn mu- 
tation in dihydropteroate synthase gene, a known target 
of thedrug (data notshown). It would be interesting if we 
could also usethe malaria mutatortodemonstrateemer- 
gence of vaccine-resistant parasites similar to the 
mutator phenotype of the bacterial pathogen, Neisseria 
meningitidis, which rapidly became resistant to bacteri- 
cidal antibody B5 and had a competitive advantage 
over a non-mutator strain. 42 
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Constructing the mutant library including various 
mutants will be useful to readily acquire novel pheno- 
types/traits. The analysed PbMut lines, however, 
shared a large proportion of variants, suggesting that 
the diverse population was not constructed after long- 
term passages in thisstudy. Mouse passage by inoculat- 
ing a large number of mutator parasites (e.g. 1 0 6 
infected erythrocytes) and increasing passage lineages 
will serve to construct the mutant library. Ifthisattempt 
improved the diversity of mutant library, more precise 
mutation hot spots might be identified by analysing 
the several number of mutant parasite lines. 

In the malaria mutator system, there isone limitation 
to overcome; unrelated mutations may be inevitably 
accumulated in the parasite during the propagation of 
parasites through mouse passages. Once the mutator 
parasite with a novel phenotype is isolated, accumula- 
tion of further unrelated mutations in the parasite 
must be stopped. Genetic crossing with the wild-type 
parasites can be used to replace the mutated DNA poly- 
merase 8 with the wild-type enzyme. Also, genetic 
crossing will be useful to reduce extra mutations unin- 
volved in the phenotypes. Although PbMut lines at the 
1 22nd mouse passage lost the ability of gametocyte 
formation, accumulation of deleterious mutations in 
genes essential for mosquito or liver stages will be sup- 
pressed by conducting a long-term mouse passage in 
combination with mosquito passage in the future. The 
candidate mutation(s) involved in novel phenotypes 
could be inferred with whole genome sequencing 
and/or linkage group selection. 6 Furthermore, screen- 
ing test with recent molecular technology like an artifi- 
cial chromosome library could be useful to accelerate 
the identification of genes involved in novel phenotypes 
of interest. 7,8 We believe thatthe mutatorsystem could 
be useful for rapid generation of novel mutants, as 
rodent malaria parasites are ideally suited for ease 
with which the whole life cycle of Plasmodium asexual 
stage could be achieved. The characterization of the 
malaria mutator will also allow rapid correlation of se- 
quence data with biological function. 
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